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JL . We find exact analytical expressions for mixing angles in matter in the con- 

M-t. 

D . 

. text of three generation neutrino oscillations in matter to discuss the role of 



X 



resonances in this phenomenon. We show that some knowledge from conven- 
tional two neutrino MSW effect, which has been extended to approximated 
solutions to three neutrino oscillations, has to be abandoned in this exact ap- 
proach. We observe that maximal values for the mixing angles in matter are 
found in nonresonant regions and stationary phases do not coincide anymore 
with resonances in this simple extension of the MSW effect. We present a gen- 
eral way to identify a resonance and discuss what we can physically expect in 
these regions. 
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I. INTRODUCTION 



Resonant regions are previleged zones for neutrino conversion. Concerning solar neutri- 
nos, the importance of a resonance can be appreciated remenbering that the standard MSW 
solution to the solar neutrino problem requires values for the mixing angle in vacuum 9 and 
for the squared mass difference A = m| — m\ such that sin 2 29 < 10~ 3 and A sin 2 29 »s 1CT 8 
eV 2 [Q which imply a resonance in the neutrino trajectory inside the sun when the approx- 
imately exponentially decreasing standard solar matter distribution is assumed ||. This is 
the so-called nonadiabatic solution to the solar neutrino problem and the role of the reso- 
nance is evident in such situation once that it is well known that the adiabaticity parameter 
@ presents its smallest values in a resonance region, which imply that neutrino transitions 
are less adiabatic in that region. 

Resonances in two family MSW effect ]3]|| are associated with maximum mixing between 
the two flavor eigenstates. This can be appreciated investigating the behavior of the matter 
mixing angle when the relevant matter density varies along the neutrino trajectory. The 
mixing angle in matter 9 is introduced as the parameter that characterizes a rotation of the 
two-dimensional neutrino space from the basis of the current eigenstates (z/ e , z/ /t ) to the basis 
of the physical eigenstates (z/ 1; i/ 2 ): 

v^t) = u e (t) cosfl(t) - i/ M (t) sin 9{t), 
v 2 (t) = u e (t) sin0(t) + 2/ M (t) cos0(t). 

It can be calculated [§]: 

sin 2 29(t) = Sin22g 2 , (2) 

' 2EV ~ 2G ^ - cos2 9] +sin 2 2fl 

where E is the neutrino energy and GpN e (t) is the consequence of electron neutrino coherent 
forward scattering from electrons in matter, the number density of which at the region 
reached by neutrinos at instant t is N e {t). 

From Eq. (^) it is possible to see that 9 is substantially modified by the neutrino coherent 
scattering from the medium. If N e (t) -^0,9^9 and we recover vacuum expressions. When 



N e (t) is extremely large, 9 — > 7r/2 and v\ — > —v^ while z/ 2 — > v e . An interesting intermediate 
case occurs when 

*-® = (3) 

and the brackets in the denominator of Eq. (|2|) vanishes. In this point the mixing of flavor 
eigenstates is maximal, i.e., from Eq. (|l|) we see that the probability of finding an electron 
or a muon neutrino in any of the mass eigenstates is 1/2. This feature has been used to 
characterize a resonance: the maximum of the bell-shaped sin 2 29 x N e graph indicates a 
resonance. 

The resonance condition given by Eq. (Q) coincides also with the position where the 
difference of the two squared matter eigenvalues of the corresponding time evolution matrix 
in matter fh\—fh\ is a minimum, suggesting that the resonance is the region where transitions 
between matter eigenstates are most likely to happen. 

Finally, it was noticed in reference |J that the resonance condition (^[) coincides also 
with the condition of existence of a stationary phase ]7]] in the two neutrino time evolution 
equations. Such fact allows to investigate the evolution of this neutrino system around a 
resonance calculating, through the stationary phase method , the related Green function. 
Employing this method it was possible to evaluate || the level crossing probability, i.e., the 
probability of nonadiabatic transitions between matter eigenstates v\ and z/ 2 as an alternative 
approach to Landau-Zener || or Petcov [[H]] methods. 

In this paper we investigate how is the behavior of mixing angles in matter and how 
to identify a resonance in the context of a three neutrino system oscillating in matter. We 
assume standard electroweak interactions of neutrinos with matter as well as nonvanishing 
vacuum mixing angles and nondegenerated mass eigenstates (in vacuum). Therefore we are 
analysing the simplest extension of the conventional MSW effect to the case where 
three families are present. We verify that the above mentioned three criteria usually used 
to define a resonance in two neutrino matter oscillations, namely, maximal mixing angles 
in matter, minimal eigenvalue difference and the presence of a stationary phase, do not 



lead anymore to the same region in the neutrino trajectory. Note also that these same 
criteria have been used in approximated solutions to three neutrino oscillations in matter 



TT| , |T2|| . Consequently some of them have to be abandoned. We present, therefore, based 
on exact analytical expressions for mixing angles in matter, how we can use our previous 
knowledge coming from two neutrino matter oscillations to arrive to a solid condition defining 
resonances in three neutrino oscillations and, therefore, an accurate analytical description 
of the physical consequences around such regions. 

II. ANALYTICAL SOLUTION 

A general time evolution equation describing a three level system can be written as an 
equation for a three-component spinor = $ 2 , ^3)'- 

i±*{t) = h(t)*(t), (4) 

where the hamiltonian h(t) is a 3 x 3 matrix which elements are specified according to the 
dynamical situation from which a boundary condition $(t ) is given. A general solution of 
Eq. (H) can be written in the the form 



= Exp 



i / h(t')dt' 



to 



HQ, (5) 



where the symbol Exp represents a sum of multiple time ordered integrals [I3|. 

For a time-independent hamiltonian, the solution of Eq. (|J) can be obtained by means 
of the Laplace transformation. Introducing the Laplace transformed ^(p) = L [$(£)], then 

p^(p) - $(t ) = -ih^ip) (6) 

and 

$(t) = L- 1 [(pl + ^)- 1 ]$(t ). (7) 

The solution Q(t) depends on the elements of the h matrix and on the roots Aj (i = 1,2,3) 
of the characteristic polynomial of the h matrix 



det [pi + ih] = 0. (8) 

In the particular case we are interested in, where a three neutrino system oscillates in 
matter, interacting with it through standard electroweak interactions, the h matrix is given 
by 

1 



h = — 
2E 



UM 2 U~ l + A\ , (9) 
where M l is a diagonal matrix given by 

(M% = rrgSa, (10) 
m\ are the three neutrino squared mass eigenvalues in vacuum, 

U = e #A Te ^A 5( >A 2 ( U ) 

is the 3x3 mixing matrix where Aj are the Gell-Mann matrices, ip, <fr and u are the mixing 
angles in vacuum and T is a matrix containing complex phases that we will ignore since we 
assume CP conservation (r = 1). 

Since we consider here only standard neutrino interactions with ordinary matter, A 
matrix has its first element An given by 

A u = 2\f2G F N e E (12) 

and all others are zero. Note that neutral current contributions to A are proportional to 
the unit matrix, giving only irrelevant overall phases to the final solution of Eq. (^J). Gp, E 
and iV e were previously introduced. 

For neutrino propagating in vacuum, A = 0, and the solution of Eq. (^) is trivial and 
simply given by 

$(t) = UrrfU-^iQ, (13) 
where m 2 is a diagonal matrix with elements 
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(m 2 ) 



'j 



exp 
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(14) 



The solution of Eq. in matter, with A being a time-dependent matrix, is given by Eq. 
(||) and it depends on the specific N e function describing the electron density. However, when 
A can be considered a constant matrix, as it is supposed in the adiabatic approximation, 
Eq. ([D has an exact analytical solution, obtained by Laplace transformation. Furthermore, 
the A matrix is invariant under a e^ A7 rotation, then, introducing now 



(t) = e"^ A7 $(t), 

we observe that ^(t) satisfies the following differential equation 

d 



-¥(t) = -iH*(t), 



with boundary condition *$f(t 



-iipA-, 



$(t ) and 



(15) 



(16) 



H 



2E 



which can be explicitly written as 



H 



AE 



A cos 2 6 + 2m\ sin 2 6 + 2A A sin 2u> cos < 



(777,3 — f ) sm 20 



A sin 2u; cos < 



V 



[m\ — |) sin! 



S + Acos2o; — Asin2u;sin0 
- -y /-Mu_< / — Asin2u;sin0 A sin 2 <fi + 2m\ cos 2 1 

where A = m\ — m 2 , Ai = m\ + m\ — 2m 2 , E = m\ + m 2 and A = S — A cos 2a;. 
On the Laplace space we have 

*(p) = [pi + iizr 1 ^). 



(17) 



;is) 



(19) 



To calculate \l/(t) we have to obtain the roots of the characteristic polynomial of the 
matrix H, det [pi + iH] = 0, which are given by |T5] 



mi + mi + mi + A 1 —Q a 

Ai = — W — cos-, 

3 3 



(20) 



A 2 



m? + r?7, + mi + A 1 / — Q a 1 

— 1 \ — -cos 

6E 2E V 3 3 2E 



-Q sin 



a 



(21) 



m? + ml + ml + A 1 /— Q a 1 / — — . a , nnS 

A 3 = — - + — \— -cos- + — J -Q sin-, 22 

3 6E 2E\ 3 3 2E y v 3 V 7 



where 



and 



a = arccos — 7==, (23) 
-1 [A 2 A? A 2 AAcos 2 0cos2cj A^cos^ - 2 sin 2 0) 1 

g= (2^F T + iT + T 2 + ^ "j (24) 



= 1 f A? _ 3 _ 9A 2 A! 3A 2 A(cos 2 0-2sin 2 0) ] 
' 27i 2 /•;)•"• [ 4 4 4 J 

1 r 9 9 

1 -A 2 A cos 2 COS 2cj - -A A 2 (cos 2 0-2 sin 2 d>) J. + 



27(2£) 3 12 r 4' 

- 27( 2 E) 3 {l AlAA cos2 4>cos2lj - Z -A 2 A x (cos 2 - 2 sin 2 0)} . (25) 
Note now that in vacuum we have 

*I = =g, = andA^g, (26) 

and the Laplace anti-transformation of Eq. (|19D reproduces the corresponding solution given 
by Eq. (|13D. The roots \ of the characteristic polynomial are the squared mass eigenvalues 
in matter. Because of the arbitrariness in the choice of the order of the roots, we use the 
above vacuum limit to order the roots in terms of the squared mass eigenvalues in the matter. 
We define: 

~ 2 ~ 2 ~ 2 

mi , mi , , m, 

Al = 2l' A! = 55 aDd = 21- (2?) 
Finally, we can write the solution of Eq. @ in terms of a T transition matrix such that 

$(£) = e^ A7 Te-^ A7 $(t ), (28) 

where the elements of the T matrix, given in terms of the Aj roots and of the elements of 
the H matrix, can be written as: 



i) diagonal elements: 



Tu — ^2 C m (A m — Hjj)(\ m — Hkk) — Hjj. 



m=l 



-iXmt 



non diagonal elements (Ty = Tji): 



(29) 



— 2_j Cm [Hij(X m — Hkk) — HikHjk] e 



-i\ m t 



m=l 



where 



(30) 



C m — [(A m — A^)(A m — A„)] 



(31) 



with m t n and n, £, n = (1, 2, 3). 

Note also that all well known results for a two neutrino system oscillating in matter can 
be straightforwardly obtained from the solution given by Eq. (p8|) . 



III. MIXING ANGLES IN MATTER 

It is well known that the knowledge of the mixing angles in matter is important to study 
resonant transitions between flavor neutrino states [|J. In order to explicitly write an exact 
expression for these angles, we define ip, 4> and uj as the mixing angles in the matter. We can 
write therefore the final solution of Eq. (f|) in terms of mixing angles in matter in analogy 
with what we did in the vacuum case, Eq.(^), using now the final solution given by Eq. 
(P8|). This solution can be written in the following way 

$(t) = U(4>,4>,Cu)M m U- 1 (4>,4>,Cu)<l>(t ) (32) 

where 

(M^ij = exp[—it\i]8ij (33) 

and 

U$,4>,£u) = e^ A V Ms e^ A2 . (34) 
8 



In order to get the matter mixing angles we simply compare Eq. (^) with Eq. (|28|) , and 
after some algebra, we obtain 



2 A x — (H 2 2 + #33) Ai + H 22 H 33 — # 23 



Sm0 " (Ai - A 2 )(Ai — A 3 ) ' (35) 



2 - [^i _ (^22 + H 33 )X 3 + H 22 H 33 - #f 3 ](Ai - A 2 ) 



tan u) v ^ 1 "ooyo i "^*-oo "-2/ (36) 

[A 2 2 — (#22 + H 33 )X 2 + H 22 H 33 — H 23 ](X 3 — Ai) 



and 



7 _ (#12 Ai + H 13 H 23 - H 12 H 33 ) cosip + (H 13 X 1 + H 12 H 23 - H 22 H 13 ) sin^ 

an ^~ (# 13 Ai + #i 2 #23 - ^22^13) COS ^ - (#12 Ai + #13#23 ~ ^12^ 33 ) sinf 



IV. ANALYSIS OF THE RESULTS 

In Fig. 1 it is presented a comparison of the behavior of the quadratic matter eigenvalues 
mf and the relevant matter mixing angles u and as a function of the parameter A for 
specially chosen values of vacuum parameters (see the corresponding caption for details. 
A is given in units of mf). There are two resonances clearly indicated by the minimum 
difference between the shown quadratic masses. We observe that u presents a maximum 
value in the lower resonance (where m| — rh\ is minimum) while shows a maximum in 
the region of the higher resonance (rh 3 — rnj is a minimum). Interesting enough, differently 
from what is expected in the two flavor neutrino oscillations in matter, the conventional 
MSW effect, and also from what was found in previous approximated analyses of the three 



neutrino oscillations JTT],|T2[, a second peak for the mixing angle u is found after the higher 
resonance [0. In Fig. 2 we show the same graphs presented in Fig. 1 to evidenced this 
unexpected behavior of the mixing angle u for larger values of A. It is clear from this figure 
that the criterion of defining a resonance by means of localizing the maximal mixing angle 
in matter, which can be safely used in two neutrino conventional MSW effect, leads to some 
ambiguity in the context of three neutrino oscillations and therefore has to be abandoned. 

Instead, we can improve this criterion analysing the content of Figs. 3 and 4. Note that 
the admixture of flavor eigenstates in each of the matter eigenstates can be obtained through 
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v i — J2 a UiaV a , where i = 1,2,3, a = e, fi,r and Ui a is given by Eq. (|TT|). Let us write now, 
as an example, the linear combination of flavor eigenstates in the first matter eigenstates: 

V\ = cos <j> cos Qv e + cos (j> sin Qu^ + sin <fiv T . (38) 

In Figs. 3 and 4 we show therefore the coefficients of this admixture (values of the vacuum 
parameters are shown in the corresponding captions). From Fig. 3 we observe that in the 
lower resonance the mixing of electronic and muonic flavor eigenstates is maximal (when 
cos (J) cos to = cos sin a)), while, from Fig. 4, we see that the higher resonance coincides with 
the maximum admixture of and v T , when cos sin d) = sin0. 

Therefore, although we detected maxima of the mixing angles in matter in regions far 
from resonances, it is still possible to identify a resonance region searching for maximal 
mixing between flavor eigenstates. Note also that such maximum are not anymore related 
with values of \/2/2 for flavor coefficients \Ui a \ in the way it happened in the conventional 
MSW phenomenon. This is because there could be nonnegligible contributions from the 
flavor eigenstate that does not participate in the resonant process. From the unitarity of 
the mixing matrix, we know that Ui e \ 2 + {Ui^ 2 + \Ui T \ 2 = 1, for i = 1,2, 3. Therefore, in the 
case where i — 1 and \Ui T \ is not vanishing, the maximal mixing is such that \U\ e \ — < 
a/2/2. A similar situation occurs for the higher resonance where we obtain \U\^\ = \Ui T \ < 
y/2/2. We can say that in three neutrino oscillation phenomenon the mixing between flavor 
eigenstate around a resonant region is as maximal as possible, although not in the same 
way as in two neutrino oscillations, where maximum mixing implies that each one of the 
neutrino flavor eigenstate participating in the resonant process contributes with 50% to the 
matter eigenstates. 

A final issue to be discussed is the criterion of identifying a resonance looking for a 
stationary phase in the neutrino evolution equations (^), in the same way it was proposed 
in reference |J in the context of two neutrino MSW effect. A stationary phase is given by 
the smallest difference of any two diagonal elements of the relevant evolution matrix when 
one of these elements is time dependent. As an example, we quote solar neutrinos where 
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the matter density considerably varies along the neutrino trajectory from the center of the 
sun, where neutrinos are created, to the solar surface. Although in two neutrino oscillations 
this criterion can be safely used, it does not work anymore in the presently analysed three 
neutrino MSW effect. Stationary phases do not coincide with the minimum squared mass 
differences or maximum flavor admixture. 

Note however that it is still possible to use the stationary phase method to calculate level 
crossing probabilities in the three neutrino oscillations. Making convenient SU(3) transfor- 
mations on the evolution matrix (|9|) it is possible to conciliate resonances and stationary 
phases. This is because resonances are invariant under similarity transformations, while 
stationary phases do not. Therefore the matrix 

H 1 = e-^He**** (39) 

presents a stationary phase for the minimum of H\ x — H\ 2 coinciding with the minimum 
of the squared mass difference rn| — rhf, and it can be used to calculate the level crossing 
probability || around the lower resonance. To obtain the correct stationary phase to analyse 
the higher resonance, we rotate the evolution matrix given in Eq. (f|) in the following way: 

H = e-^ A7 he^ A7 (40) 

Now the minimum of the difference H 33 — Hu indicates a stationary phase which now 
coincides with the required resonance. 

V. CONCLUSIONS 

Resonances represent a crucial region in the time evolution of neutrinos oscillating in 
matter. They are closely related with the nonadiabatic character of the oscillation. We 
investigated a general criterion to define a resonant region when three neutrino are present 
in the oscillation phenomenon. We observed that two of the three commonly employed 
criteria to identify a resonance in two neutrino oscillations are not valid anymore in its 
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simplest extension to three neutrino MSW effect. For instance, mixing angles can present 
maximal values far from resonant regions and therefore this criterion to define a resonance 
has to be abandoned. Furthermore, stationary phases do not necessarily coincides with 
resonant regions. The safest way to identify such resonance regions is to investigate the 
behavior of the squared matter eigenvalue differences, looking for their minimum values. 



VI. ACKNOWLEDGEMENTS 

The authors would like to thank Fundagao de Amparo a Pesquisa do Estado de Sao 
Paulo (FAPESP), Conselho Nacional de Desenvolvimento Cientifico e Tecnologico (CNPq) 
and Coordenagao de Aperfeigoamento de Pessoal de Nivel Superior (CAPES) for several 
financial supports. 



12 



REFERENCES 

[1] V. Barger, R.J.N. Phillips and K. Whisnant, Phys. Rev. D43 (1991) 1110; 

[2] J.N. Bahcall and R. Ulrich, Rev. Modern Physics 60, 297 (1988). 

[3] S.P. Mikheyev and A. Yu. Smirnov, Sov. J. Nucl. Phys. 42, 1913 (1985); Nuovo Cimento 
9C, 17 (1986). 

[4] S.M. Bilenky and S.T. Petcov, Rev. Mod. Phys. 59, 671 (1987). 

[5] L. Wolfenstein, Phys. Rev. D 17, 2369 (1978); 20, 2634 (1979). 

[6] M.M. Guzzo and J. Bellandi, Phys. Lett. B 294, 243 (1992). 

[7] A. Erdelyi, Asymptotic Expansions, Dover Pub. Inc. NY (1956). 

[8] M.M. Guzzo, J. Bellandi, and V.M. Aquino, Phys. D 49, 1404 (1994). 

[9] C. Zener, Proc. R. Soc. London A 137, 696 (1932); L. Landau, Phys. Z. Sow. 2, 46 
(1932). 

[10] S.T. Petcov, Phys. Lett. B 191, 299 (1987). 

[11] T.K. Kuo and James Pantaleone, Phys. Rev. D 35, 3432 (1987). 

[12] V. Barger, K. Whisnant, S. Pakvasa, and R.J.N. Phillips, Phys. Rev. D 22, 2718 (1980). 

[13] Some physical consequences of this fact have been analysed in M.M. Guzzo, O.L.G. 
Peres, V. Pleitez and R. Zukanovich Funchal, "Nonresonant MSW solution to the solar 
neutrino problem", preprint IFT-P.003/95, IFUSP-P/1133. 

[14] R.P Feynman, Phys. Rev. 84, 108 (1951). 

[15] Harry Bateman, Tables of Integral Transforms, Mc Graw-Hill, Book Company Inc. 
(1954). 



13 



CAPTIONS 



Figure 1: Squared matter eingenvalues m 2 and the relevant matter mixing angles Co and 
as a function of the parameter A are presented. The values of the parameters 0, oj, and m 2 in 
vacuum: m\ = 5m 2 , = 25m 2 ; sin 2 = 5 x 1CT 4 ; sin 2 oo cos 2 = 5 x 1CT 2 were chosen in order 
to well demonstrate the behavior of these parameters as a function of A. The eingenvalues 
rhi and the quantity A are given in units of m 2 . 

Figure 2: Mixing angles Co and as a function of the A for larger values of A. The 
parameters (f>,iv, and m 2 are the same as that of the Figure 1. 

Figure 3: The squared mass difference m 2 — m 2 , and the quantities cos sin Co and 
cos cos Co quantities are presented as a function of the energy of neutrinos. The values 
of the parameters (f>,u> and mf in vacuum are: m| = 1.445 x 10~ 4 eK 2 ; m\ = 10~ 8 eV 2 ; 
m\ = 0; sin 2 = 5 x 10~ 4 ; sin 2 oj = 0.050025. The squared masses are given in units of m\. 

Figure 4: The squared mass difference m\—m\ , the effective cos sin Co and sin0 quantities 
are presented as a function of the energy of neutrinos. The values of the parameters , uj 
and mf in vacuum are the same ones used to draw Fig. 3. 
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